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^ ■ Abstract 

The diquark model is used to calculate the electromagnetic polar- 
izabilities and charge radii of the nucleons for three different poten- 
tials. Making the scalar diquark lower in mass introduces a mixing 
angle 9 between the |56) and |70) states ,which allows an improve- 
ment in value of all 6 properties. Generalizing the Gamov- Teller ma- 
trix and the magnetic moment operator to the diquark model gives 
constraints on this mixing. We obtain for the Richardson potential 
e = 23.2°, ap = 7.9lJ;0 x IQ-^fn^^ a„ = 7.7l0j x IQ-^fm^ ^^ = 
5.4+Ji X 10-Vm3, A = G-7tl-'r x 10-^^^, (r^)^ = 0.37t^:°^/m^ 
(r^) = — 0.07lQ'Q2/m^. Additional pion cloud contributions could 
improve on all six results. 
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The electric and magnetic polarizabilities, labeled a and (3 , respectively, 
have been measured recently [j^]'[^, yielding the results Op = (10.9 ± 2.2 it 1.3); 
10-^ fm^, 'Pp = (3.3=f2.2=f1-3) x 10-V"^^ «„ = (12.3 ± 1.5 ± 2.0) x 
10~^/m^ and /?„ = (3.5 =F 1.5 =F 2.0) x 10~''/m^. The experimental values 
are obtained by fitting the Compton scattering data to relations obtained 
from low energy theorems Q. In addition to the polarizabilities we will 
be concerned here with the charge radii (r^) = ^Cixf, whose values were 
measured to be|] (r^) = O.TOS/m^ and (r^)^ = -0.113 ± 0.003/m2. The- 
oretical studies, using the 3 symmetric quark modeip, Skyrme modelsjg], 
MIT-bag models P] or nonlinear meson theories P] have not obtained satisfac- 
tory results for all these quantities. A recent paper using chiral perturbation 
theory]^ gives good agreement for a, but predicts negative values for the 
magnetic polarizabilities. 

In this paper we will use another approach, the diquark model. This 



model, which provides improved agreement with a wide range of data|lC], 
reduces the three body problem of the quark picture to a two body problem 
and results, therefore, in a considerable simplification of the calculation of 
baryonic properties. To review the salient features, the proton and neutron 



wavefunctions have spin-flavour forms 
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wliere S and t are two different diquark states (S for sextet and t for triplet 
of SU(3) flavour) with spin S=l and S=0, respectively. The parameter a=l 
for the fully SU(6) symmetric scheme. The spin interaction between the 
two quarks, which form the diquark, is assumed to yield a mass splitting 
between the S- and the t-state and breaks SU(6) symmetry. From the QCD 
hyperfine interaction Am (xSi ■ S2 it follows that the S-state is heavier than 
the t-diquark. Hence a will differ from 1, as we discuss below. In addition 
there will be another mass splitting A?n,2 coming from the spin interaction 
between the quark and the diquark, which increases the spin | mass relative 
to the spin ^ state. 

To determine the spatial dependence of the wavefunctions we use three 
different forms of scalar potentials, a logarithmic one Vl (^) = C • Log ( ^ ) , 
a superpostion of linear and Coulomb-potential Vcl (f) = — ^ + C^^ and 



the Richardson potential ||ll[] Vr, which interpolates the Coulombic potential 



for small r and the linear potential for large r. 

The expressions for the polarizabilities, derived for the nuclei of atoms 0] 
can be modified for the nucleons, yielding 

(y = —^O^e,x,) + — ^ ^ (3) 
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where a is the fine structure constant, M the mass of the nucleon, mj the 
mass of its constituents, D the electric dipole operator and /Uj the magnetic 
moment operator. 

The nonrelativistic wavefunctions for all three potentials were solved 
numerically. To get the parameters of the different potentials we fitted the 
eigenenergies to different baryon states, the N(938), N(1710), A(1232) and 
A(1600). In the diquark picture [l^] they correspond to the (56,0), (56,0)** 
(56,0) and (56,0)*, respectively, where the A wavefunction 
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has to be used for the two latter states. The dipole operator requires that the 
intermediate states of the electric polarizabilities eq.(y) are all L=l states. 
Therefore we have to take (56,1) and (70,1) states as well as their excitations. 
The 1 56) states have the same form as eqs.(||) @ and the |70) states are 



proportional to|13] 
|70+)= l^Suud-^Sudu] (vf|l'i)-VilO'T)j - \lltudu\0,])\ 

Excited states soon give a negligible contribution to a, since the denom- 
inators En — Eq increase as the numerator decreases. It is sufficient to take 
the 2p and 3p wavefunctions into account. As we will see below it is nec- 
essary to introduce a mixing between the |56) and the |70) states. To see 
the essential features of a "symmetric diquark model" in what follows we 
would treat the nucleons as if they were pure |56) states (i.e. set a = 1 in 
the wavefunctions). Defining Cu,s = m-Zn s ^'^'^ '-^2i,s = rn+mt s ^^ ^et for 
the electric polarizabilities 
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where r^^j are the coordinates of the reduced mass. 

For the magnetic polarizabihties the first term of eq. (13) needs some 
modification. The magnetic moment operator has to be generahzed for the 
quark-diquark states to 

/., = g^S^ + ^a^ (7) 

However this alone leads to problems when calculating the ratio of the mag- 
netic moments — . Only a value of g=0 would reproduce the successful 
predictions of the symmetric quark model of —1.5 (compared to the ex- 
perimental value of —1.46). Analyzing the difference between the diquark 
picture and the latter reveals the need for a contribution coming from 
{Sud\ fJ-z \tud) , which is non-zero for a composite diquark. Therefore we add 
{Sud\ 2^0-"^ + ■^cr'' \tud) to eq.(0). If we include this term we get exactly 
the quark result again, in the limit g— > 2 and rrid — > 2m. This shows how 
important it is not to see the diquark naively as a single point particle but 
rather as a composite of two quarks. To continue, it is sufficient to take the 



A as the dominant intermediate state [14] in eq.(Q) and neglect the higher 
mass states. The second and third terms of eq.(^ are straightforward to 
calculate and so the magnetic polarizabihties become 
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The charge radii in the diquark picture are 



vt) (r^2 I 0/^2 \ I Vs) 



(9) 



(r^) = ^(^ («^ (-?) (C'f. - Cl) + (r,^) (CL - Cl)) (11) 

Using these formulas with a = 1 we obtain the results shown in table 

1. The values of the electric polarizabilities are too small, the magnetic 

quantities too big and the charge radii too small. Note especially that the 

neutron charge radius has the wrong sign and is very close to zero. 

So far we have treated the nucleons as pure |56) states. The situation 
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Logarithmic 


Lin. -Coulomb. 


Richardson 


ap [fm^] xlO~^ 


4.1 ±0.4 


5.0 ±0.4 


5.2 ±0.5 


an [fm^] xlO-4 


3.3 ±0.1 


4.5 ±0.1 


4.4 ±0.2 


:5p[/m3] xlO-4 


3.9 ±0.2 


7.8 ±0.4 


7.6 ± 0.2 


:5„ [/m3] xlO-4 


4.1 ±0.2 


8.4 ±0.4 


9.0 ±0.4 


(-')p [/-'] 


0.20 ±0.02 


0.25 ±0.01 


0.25 ±0.02 


(-^)„ [/-^] 


0.005 ± 0.004 


0.012 ±0.006 


0.015 ±0.006 



Table 1: Results for the polarizabihties and charge radii for different forms 
of potentials with no mixing, a=l 

turns out quite favourable if we introduce a mixing angle 9 between the |56) 
and the |70) states to form the nucleons 



\N) = cos 61 1 56) + sin 61 1 70) 



(12) 



which is equivalent to introducing a mixing parameter a, for which cos 9 = 
,"^ . The QCD hyperfine splitting as well as phenomenological studies [p^] 
suggest that the S=0 state has lower mass and thus is a larger component 
of the nucleon. We will use values a > 1 in the above equations. To get an 
upper constraint on a, we fit our model to the ratio of the nucleon magnetic 
moments — as well as to the axial to vector current coupling constant ratio 

gv ' 



Using the above magnetic moment operator for the diquark model and 

the wavefunctions eqs.(|^),(^) we obtain 

f a'^ + a g \ 

— = -3^4 —^ 13 

\ 2m rrig I 

The operator for ^, also called the Gamov- Teller matrix|17] is 

(p,TlEir+5f|n,T) 
(p,TlE.r+|n,T) ^ ' 

As we generalized the magnetic moment operator to the diquark model we 
will also have to generalize the Gamov- Teller matrix. Analogously we start 
with 

i 

h can be factorized as h = t*g and comparing the terms to the quark picture 
shows that t= -y=. As in the case of the magnetic moments we have to add 
terms from the quark picture, 

(r+cjf + T+al) |tw, 0) = V2\Suu, 0) 

(tiX + r2+a|) \Sdd,0) = V2\tud,0) 

The denominator of eq.([lj) remains —1 also for the case of mixing and 
so we get 

f;^-(.,T,i:.^SfKT,.l±|^±^ ,15, 



Now we take the ratio of the magnetic moments fixed at its experimantal 
value — = — 1.46p^] and use eq.(13) to calculate g as a function of a and 



^, which in turn (eq.(p^)) gives ^ as another function of those parameters. 
The dependence on ^ is weak compared to the a-dependence. A value of 



a = 2.5 ± 0.1 fits the experimental value of ^ = 1.25[18]. Since taking the 



relativistic nucleon wavefunctions would have a lowering effect anyway on 



^[19| we will take a = 2.5 as the upper limit and take Aa = —0.3 for a 



theoretical error calculation. In terms of the mixing angle of eq.(|12D a = 2.5 
means 9 = 23.2° and therefore gives a small |70) state contribution to the 
nucleons. 

Figure 1 shows the results for our quantities as a function of the mixing 
factor a for the Richardson potential. The situation for the other potentials 
looks qualitatively similar. We see immediately that a large mixing factor 
improves the data for all 6 quantities and justifies our approach. 

Table 2 shows the results for all 3 forms of potentials, where the diquark 
mass splitting was taken as Am = 100 it 50MeV. Although the absolute 
values for the charge radii are still too small they acquired the right sign. 
The results for the superposition of the linear & Coulombic potential and 
the Richardson potantial lie within the error bars of the the experimen- 
tal results. The magnetic polarizabilities are on the upper limit and the 
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ap [fm^] xlO-^ 


tr O+0.4 
^•0-0.8 


8.4il:i 


7 q+1-0 
'•^-0.9 


an [/m3] xlO-4 


c: 9+0.4 


8.4+0:i 


7 7+0.3 
'•'-0.6 


:5p[/m3] xlO-4 


2.0l°i 


fi t;+1.0 
0-^-0.6 


c; 4+1.6 


:5„ [/m3] xlO-4 


2.7tL' 


7 q+2.0 
'•^-0.9 


fi 7+1-3 
"•'-0.7 


(-')p [/-'] 


0.271°:°^ 


o.4oj:°:°^ 


n 07+0.02 
'-'•'^'-0.03 


(-^)„ [/-^] 


-o.oeitnS^ 


-0.08tn2 


-n 07+0-03 

'J-'-" -0.02 



Table 2: Results for the polarizabilities and charge radii for different forms 
of potentials 

electric quantities on the lower edge. The magnitutes of the charge radii 
are too small by about 50%. Overall, the ability of the diquark model to 
give reasonable values for all the static properties of the nucleon are very 
encouraging. 

Analyzing the contributions of the different terms shows us that an in- 
crease of the terms containing the spatial wavefunction would lower the 
electric and raise the magnetic values of the polarizabilities. A pion cloud 
around the nucleon core gives exactly such a contribution. Such sources 
were already included by Weiner and Weise pq] and had considerable effects 
on the electric polarizabilities. Pion clouds will also effect the charge radii. 
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The transitions p ^ n + vr"*" give positive contribution and n ^ p + tt~ 
yields the desired negative part. Therefore we have begun to consider the 
influence the pion cloud will have on the constraints — and — and to find 
the corresponding values for the mixing a. The core values can be directly 
taken from this paper - the pion cloud terms will have to be added. 
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Figure Captions 

1. Electric (a), magnetic (b) polarizabilities and charge radii (c) for 
neutron (solid line) and proton (dashed line) as a function of the mixing 
parameter a, which is explained in the text. 
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